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Abstract 

Scintillation properties of pure Csl crystals used in the shower calorimeter being 
built for precise determination of the 7r + — >7r°e + f e decay rate are reported. Seventy- 
four individual crystals, polished and wrapped in Teflon foil, were examined in a 
multiwire drift chamber system specially designed for transmission cosmic muon 
tomography. Critical elements of the apparatus and reconstruction algorithms en- 
abling measurement of spatial detector optical nonuniformities are described. Re- 
sults are compared with a Monte Carlo simulation of the light response of an ideal 
detector. The deduced optical nonuniformity contributions to the FWHM energy 
resolution of the PIBETA Csl calorimeter for the TT + ^te + v 69.8 MeV positrons and 
the monoenergetic 70.8 MeV photons were 2.7% and 3.7%, respectively. The upper 
limit of optical nonuniformity correction to the 69.8 MeV positron low-energy tail 
between 5 MeV and 55 MeV was + 0.2 %, as opposed to the + 0.3 % tail contribution 
for the photon of the equivalent total energy. Imposing the 5 MeV calorimeter veto 
cut to suppress the electromagnetic losses, GEANT-evaluated positron and photon 
lineshape tail fractions summed over all above-threshold ADCs were found to be 
2.36±0.05(stat)±0.20(sys)% and 4.68± 0.07(stat)± 0.20(sys) %, respectively. 
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1 Introduction 



The PIBETA collaboration has proposed an experimental program [1] with 
the aim of making a precise determination of the 7r + ^7r°e + z/ e decay rate at the 
Paul Scherrer Institute (PSI). The proposed technique is designed to achieve 
an overall level of accuracy in the range of ~0.5%, improving thus on the 
present branching ratio uncertainty of 4% [2]. The (1.025± 0.034) xl0~ 8 pion 
beta decay (nfi) branching ratio will be remeasured relative to the 10 4 times 
more probable -n + ^>e + u e decay rate, that is known with the combined statis- 
tical and systematic uncertainty of ~0.40% [3,4]. 

The Standard Model description of the lift decay and its radiative correc- 
tions [5] enables a stringent test of the conserved vector current (CVC) hy- 
pothesis and the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix. The most accurate extraction of the CKM matrix element V u d, 
based on the superallowed Fermi transitions in nuclei, involves the theoretical 
nuclear overlap corrections that dominate the 0.3% total uncertainty [6]. Re- 
cent measurements of nine different nuclear transition rates confirm the CVC 
hypothesis at the level of 4xl0~ 4 , but violate the three-generation CKM uni- 
tarity condition by more than twice the estimated error [7]. Complementary 
neutron decay experiments yield a result that also differs from unitarity by 
over two standard deviations, but in the opposite sense [7]. It is of fundamen- 
tal interest to measure the branching ratio of the 7r + ^7r°e + z/ e decay at the 
0.5 % level or better to check the consistency with the nuclear and neutron 
j3 decay values and the unitary prediction of the Minimal Standard Model 
(MSM). Moreover, the precisely measured 7i[3 branching ratio could be used 
to constrain masses and couplings of additional neutral gauge bosons in grand 
unified theories. If such couplings are directly discovered at future collider 
experiments, precision low-energy data will be an essential ingredient in ex- 
tending the MSM framework [8]. 

The central part of the PIBETA detector system is a high-resolution, highly 
segmented fast shower calorimeter surrounding the active stopping target in a 
near-spherical geometry. The requirements imposed on the calorimeter design 
were: 

• high energy resolution for effective suppression of background processes; 

• high segmentation and fast timing response to handle the high event rates; 

• acceptable radiation, mechanical and chemical resistance; 

• compact geometry to simplify the operation and reduce the cost. 

These design constraints are met in a spherical detector with individual calo- 
rimeter modules made from undoped cesium iodide (Csl) crystals. The nonac- 
tivated alkali iodides have been known for more than thirty years to exhibit a 
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near-ultraviolet emission component when excited by ionizing radiation [9,10]. 
The pure Csl material was reintroduced as a fast, rugged, and relatively inex- 
pensive scintillator material by Kobayashi et al. [11] and Kubota et at. [12,13]. 
The scintillation characteristics of pure Csl crystals have been reported subse- 
quently by a number of experimental groups. These investigations have covered 
light readout techniques, scintillation decay times, origins of different emission 
components, crystal light yields, and energy and timing resolutions [14-19], 
radiation resistances [20,21], wrapping and tuning methods and the unifor- 
mity of light responses [22,23]. The design and performance of pure Csl elec- 
tromagnetic calorimeters have been recently described in Refs. [24] (NMS at 
LAMPF/BNL), [25,26] (KTeV at Fermilab), and [27] (PHENIX at RHIC). 

The PIBETA calorimeter geometry shown in Fig. 1 is obtained by the class 
II geodesic triangulation of an icosahedron [28]. Selected geodesic breakdown 
results in 220 truncated hexagonal, pentagonal and trapezial pyramids cov- 
ering the total solid angle of 0.77 x An sr. Additional 20 crystals cover two 
detector hemispheres open to the beam and act as shower leakage vetoes. 
The inner radius of the crystal ball is 26 cm, and the axial module length is 
22cm, corresponding to 12 Csl radiation lengths (X =1.85cm [29]). There 
are nine different module shapes: four irregular hexagonal truncated pyramids 
(we label them HEX-A, HEX-B, HEX-C, and HEX-D), one regular pentag- 
onal (PENT) and two irregular half-hexagonal truncated pyramids (HEX-D1 
and HEX-D2), and two trapezohedrons which function as vetoes (VET-1, 
VET-2). The volumes of our Csl crystal detector modules vary from 797 cm 3 
(HEX-D1/2) to 1718cm 3 (HEX-C). Dimensioned drawings of the HEX-A 
and HEX-D 1 crystal shapes are shown as examples in two panels of Fig. 2. 

All key components of the complete PIBETA detector have been prototyped 
and built and have met the required specifications [30]. All the detector com- 
ponents were delivered to PSI by mid-1998. The final assembly of the appa- 
ratus was followed by the in-beam commissioning and calibration. "Produc- 
tion" measurements are scheduled for 1999, and will extend over two more 
years to complete the first phase of the project, a ~ 0.5 % measurement of the 
7r + ^7r°e + z/ decay rate. 

The quality of the delivered Csl calorimeter crystals was controlled in the 
following program of measurements: 

(1) check precisely the crystal physical dimensions against the specifications; 

(2) study the effects of different crystal surface treatments and wrapping 
configurations on the energy and timing resolution and on the uniformity 
of light response; 

(3) determine the percentage of the fast scintillation component in the light 
output; 

(4) measure the contribution of detector photoelectron statistics to energy 
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resolution; 

(5) find the temperature dependence of a detector's ADC readings; 

(6) determine axial and transverse optical nonuniformities for each crystal. 

The ultimate goal of the program was to provide the input parameters for a 
Monte Carlo simulation of the calorimeter response and to deduce the effects of 
Csl crystal light yields and optical nonuniformities on the energy lineshapes 
of detected photons and positrons ranging in energy between 10 MeV and 
120 MeV. 

The spectra generated by cosmic muons have been used to examine the unifor- 
mity of light response of two large (24x36 cm) cylindrical Nal(Tl) calorimeter 
crystals in the earlier work of Dowell et al. [31]. Their apparatus relied on 
plastic scintillator hodoscopes with 3.8x3.8 cm 2 cross section for the charged 
particle tracking. The precision of the trajectory reconstruction was therefore 
limited to ~4cm. Such positional resolution was deemed inadequate for our 
application because the PIBETA Csl crystals have an axial length of 22 cm 
with front (back) surface side lengths of less than 4 cm (7cm). 

Most of our measurements were done with the tomography apparatus built 
around three pairs of multiwire drift chambers (MWDC) using cosmic muons 
as the probe. One dozen PIBETA Csl detectors were examined in a 350 MeV/c 
minimum ionizing muon beam of the PSI 7rEl well as in the 405 MeV / c 

penetrating pion and stopping proton beams in the 7rMl channel. Selected 
Csl crystals were also scanned with a 662 keV 137 Cs 7-source. The light yield 
nonuniformities measured in the beam and with a radioactive-source scan 
method were compared with more precise cosmic muon tomography data. 
The experimental apparatus used for a 7-source crystal scans and the subse- 
quent data reduction were simple enough to be easily adopted by the crystal 
manufacturers, thus accelerating the subsequent quality control cycle. 



2 Crystal Production, Mechanical Quality Control and Surface 
Treatment 

Pure Csl crystals produced for the PIBETA calorimeter came from two dif- 
ferent sources. Twenty-five crystals were grown and cut to specified shapes in 
the Bicron Corporation facility in Newbury, Ohio. The rest of the scintilla- 
tors were grown in the Institute for Single Crystals in Kharkov (AMCRYS), 
Ukraine. Preliminary examination of fifty-five AMCRYS Csl crystals, includ- 
ing the mechanical and physical tests, was done in the High Energy Institute 
of Tbilisi University, Georgia. 

Manufacturing tolerances of the crystals were specified for the linear dimen- 
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sions (+150 lira/— 50 lira) and for the angular deviations (+0.040°/— 0.013°). 
Geometrical dimensions of the machined crystals were measured upon deliv- 
ery at PSI using the computer-controlled distance-measuring device WENZEL 
Precision. The machine was programmed to automatically probe the surfaces 
of a subject crystal with a predefined shape. Each crystal surface plane was 
scanned with a touch head at six points and the equations of planes were 
found through these surveyed points. Body vertices, measured with an abso- 
lute precision of 2 /im and reproducible within 20 /im, were then compared 
with the expected theoretical values. Those crystals that failed the imposed 
geometrical tolerances were returned to the manufacturer for reuse as raw 
crystal growing material. 

After physical measurements crystal surfaces were polished with a mixture of 
0.2 [im aluminum oxide powder and etylenglycol. Next came a measurement of 
the light output of the fast Csl scintillation component (F) that is completely 
decaying in the first 100 ns, relative to the total Csl signal (T), integrated in 
a 1 /is ADC gate. These measurements were made with unwrapped crystals 
using air-gap coupled photomultiplier tubes (PMTs) and a Tektronix TDS 
744 digital oscilloscope. Only the crystals with a fast-to-total component ratio 
(F/T) better than 0.7 were accepted. The mean value of fast-to-total ratio for 
all accepted Csl crystals was 0.788, as shown in Fig. 3. 

EMI photomultiplier tubes 9821QKB [32] with 75 mm diameter cathode were 
glued to the back faces of hexagonal and pentagonal Csl crystals using a 
300 lira layer of silicone Sylgard 184 elastomer (Dow Corning RTV silicon 
rubber plus catalyst). Smaller half- hexagonal and trapezial detector modules 
were equipped with two inch EMI 9211QKA phototubes [32]. The resulting 
crystal-photomultiplier tube couplings were strong and permanent, but could 
be broken by application of a substantial tangential force. The PMT quartz 
window transparency, peaking at ~ 380 nm [32] , is approximately matched to 
the spectral excitation of a pure Csl fast scintillation light component with 
a maximum room temperature emission at ~310 nm [15]. The PMT high 
voltage dividers were modified EMI-recommended bases designed and built 
at the University of Virginia. These dividers minimized the so-called "super" - 
linearity exhibited by many PMTs well below the onset of saturation. The 
maximum PMT nonlinearity measured with a pair of light-emitting diodes 
was less than 2 % over the full dynamic range expected in the nf3 decay rate 
measurement [33]. 

Four different wrapping materials for lateral crystal surfaces were investigated 
in wrapping and tuning studies with 405MeV/c pion and proton beams in 
PSI (iMl area: 

(1) one to five layers of a 38 lira PTFE Teflon sheet (CF2 monomer, Xo=16.0 cm); 

(2) one to five layers of 30 fira Mylar sheet (C 5 H 4 2 , X =28.7cm); 
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(3) a 250 /iin thick Tyvek fleece (polyethylene CH2, Xo=131.3cm); 

(4) a 110 fim thick Millipore filter (polyvinylidene fluoride PVDP, X =19.0 cm) 
with 0.22 fj,m pores; 

(5) a wavelength shifter lacquer treatment of crystal surfaces. 

In all cases a primary diffuse reflector was protected by an additional 20 /im 
thick aluminized Mylar cover. These preliminary tests showed more uniform 
axial light collection from the crystal front section when a polished front crys- 
tal surface was covered with a black paper sheet. A similar tuning method, 
using black paper strips at the front section of 5x6x36 cm 3 CsI(Tl) crystals 
in order to improve the uniformity of light collection, was recently described 
in Ref. [23]. 

The maximum fast-to-total ratio was always achieved with unwrapped crys- 
tals because all applied wrapping materials absorb more ultraviolet light which 
dominates the fast component, than visible light. A two-layered Teflon cover 
was found to be superior in delivering >20 % more fast scintillation light com- 
ponent than the other reflectors. It was also comparatively better in not de- 
grading the F/T ratio by more than ~2%. These differences in measured 
light output were reproducible in the beam tests with ~ 1 % event statistics 
and estimated 2 % systematic uncertainty. They are subsequently confirmed 
in more controlled cosmic muon tomography measurements. 

On the basis of these findings the adopted standard wrapping configuration 
for all studied Csl crystals consisted of: (a) the lateral surfaces being wrapped 
with two layers of Teflon foil plus one layer of aluminized Mylar sheet, and 
(b) the front crystal surface covered with the black paper template. The back 
crystal surface with the glued phototube was left uncovered. 

Light yield measurements were repeated for a subset of crystals after a six 
month period to confirm that no appreciable changes occurred as a result of 
degradation of the surface reflectivities and wrapper material qualities. The 
measured light output was typically within 5 % of the originally measured 
value. 

The most successful surface treatment of Csl crystals, however, involved paint- 
ing the lateral crystal surfaces with a Waveshift lacquer. We recently studied 
this method on a large sample of Csl detectors. The light response and unifor- 
mity properties were noticeably improved, resulting in ~ 20 % better energy 
resolutions for the 70 MeV monoenergetic e ± and 50-82 MeV tagged photons. 
The possible degradation and change in the crystal surface reflectivities that 
could be difficult to account for in the multi-year long precision experiment 
like the PIBETA is also expected to be arrested by such a treatment. These 
measurements will be covered in detail in a forthcoming paper [34]. 
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3 Tomographic Apparatus 



A simplified sketch of the tomographic apparatus layout is illustrated in a 
GEANT rendering showing a few simulated cosmic muon trajectories on Fig. 4. 
Three identical delay-line readout drift chambers were used to define the cos- 
mic muon tracks intersecting Csl crystals. The chambers were built at the 
Los Alamos Meson Physics Facility (LAMPF) by the M. Sadler group from 
the Abilene Christian University. This type of drift chambers has been used 
reliably for years with several LAMPF spectrometers [35-38]. Each cham- 
ber consists of two signal and three ground planes with the nominal active 
area of 60x60 cm 2 . The horizontally oriented u x" and "y" signal planes are 
two orthogonal sets of alternating cathode and anode wires evenly spaced at 
0.4064 cm. 

The light-tight aluminum box had two drift chamber pairs mounted on its 
top plate, and one pair fixed below the bottom side. Distances between the 
centers of chamber pairs 1-3 and 2-3 were 24 cm and 27 cm, respectively. The 
dark box could accommodate up to six Csl calorimeter modules at one time. 
It had feedthrough connectors for twelve signal and twelve high voltage cables 
as well as six temperature sensor lines and six LED calibration signal cables. 
A pair of 1 cm thick plastic scintillators were placed directly below the appa- 
ratus and separated from the chamber system by a 5 cm layer of lead bricks 
shielding out the soft cosmic ray component. The Monte Carlo simulation (see 
Sec. 8.1) showed that hard cosmic muons penetrating the frames of the appa- 
ratus, Csl crystals and shielding material and triggering two scintillators had 
a smooth energy threshold starting at ~ 120 MeV. The zenithal angle range of 
accepted cosmic muon trajectories intersecting all three MWDCs and at least 
one Csl detector was ±45°. 

The chambers were operated with anode wires held at positive voltage between 
2400-2600 V. The chamber cathode wires were grounded. The gas mixture was 
65 % argon, 35 % isobutane plus 0.1 % isopropyl alcohol. Detection efficiencies 
of individual chambers for penetrating cosmic muons exceeded 90 %, with the 
combined six-chamber efficiency routinely surpassing 50 %. 

Anode wires were attached directly to a fast 2.5 ns/cm delay-line. Signals 
from both ends of the delay line were amplified twenty-fold, discriminated 
at a threshold of 10 mV and connected to two channels of a time-to-digital 
converter (LCR 2229A TDC). The time difference in the two TDC readings 
identified the fired anode wire. Cathode lines defining the alternate field [39,40] 
were connected to one "odd" (O) and one "even" (E) line, and the signals on 
these lines were processed by a custom-made electronics unit which added 
and subtracted the analog pulses [41]. The electronic sum of cathode pulses 
(O+E) was discriminated and the resulting delayed signal determined the drift 
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position timing, and was used to stop another TDC channel. The difference of 
the cathode signals (0— E) was digitized with an analog-to-digital converter 
(LCR 2249A ADC). That information was used to discriminate between the 
events that produced the ionization tracks left and right of the given anode 
wire. The 100 ns integration gate timing was defined by the (O+E) logic pulse 
timing. Four measurements were therefore required to find an intersection 
between a cosmic muon track and one drift chamber plane, namely three TDC 
values and one ADC value. The schematic diagram of the electronic logic is 
shown in Fig. 5. 

Advantages of the system were good charged particle track resolution (~ 
0.5mm root-mean-square in the horizontal plane, Fig. 6), stability to fluc- 
tuations in the outside environmental parameters (humidity, temperature and 
pressure) and low cost of the associated electronics logic and readout. Lim- 
ited counting rate and restriction to single hit events in the chambers did not 
represent a drawback in this tomographic application. The apparatus was op- 
erated in an air-conditioned room with a controlled humidity level kept below 
30 %. The temperature at six points inside the dark box as well as absolute 
time were recorded for every triggered event. The temperature range recorded 
inside the dark box during three years of data taking was (22°±3°)C. The 
typical temperature gradient was 0.4° C/day, leading to average absolute tem- 
perature variation of 0.2° C and average Csl ADC gain drift of ~0.4% in a 
single data taking run (fixed at 250000 triggers, ~6 hours). Measured light 
output variations of the Csl crystals and the PMT gain drifts caused by daily 
temperature cycles were compensated for in the replay analysis (see Sec. 6). 



4 Data Acquisition System 



The computer code used for data acquisition was HIX (Heterogeneous Infor- 
mation Exchange), a data acquisition system developed by S. Ritt at PSI, 
originally intended for use in small and medium-size nuclear physics experi- 
ments [42]. 

The "frontend" 486 personal computer was connected to a CAMAC crate via 
a HYTEC 1331 interface that read ADC, TDC, scalar, and temperature sen- 
sor units. The C program running under MS-DOS accessed data and sent it 
over the network to a VAX 3100 server. A simple communication program on 
a VAXstation received data from the frontend computer and stored it in a 
global section buffer. Buffered data were passed by the Logger application to 
a user-written Analyzer program. Analyzed events were subject to predefined 
cuts, filled in predefined histograms, and stored in a raw data stream directed 
simultaneously to a hard disk drive and an 8-mm tape system. The exper- 
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iment was controlled from a PC computer by a Microsoft Windows control 
program. A windows-based graphical user interface to the program allowed 
starting, pausing and stopping data acquisition, as well as online inspection 
of individual raw data words, calculated data words, scalers, assorted efficien- 
cies, one and two-dimensional histograms, and gate, box and Boolean tests. 
Different PC computers on the Ethernet could make a connection to the VAX 
server called Link and access the same information remotely. 

During data acquisition only those events for which at least two x and two y 
drift chambers had nonzero ADC and TDC values were written to an disk- 
resident ASCII Data Summary File (DSF). One good event contained 55 raw 
data words specifying the response of the MWDCs and Csl detectors, instan- 
taneous temperatures and absolute times. That was a full set of observables for 
which a cosmic muon trajectory could be unambiguously reconstructed. The 
raw trigger rate was ~ 13 Hz with six Csl crystals in the dark box. The average 
DSF event rate was ~ 4 Hz. Individual runs were stopped and restarted auto- 
matically after 250,000 collected triggers. A total of about 10 5 DSF events per 
crystal were typically collected in one week of data acquisition. Tomography 
data for one set of six crystals were usually collected over a two week period. 



5 Drift Chamber Calibration 



Two different trigger configurations were used interchangeably during the data 
collection. Drift chamber calibration was done with a trigger requiring a two- 
scintillator coincidence and good ADC and TDC data for at least one drift 
chamber. In the tomography data acquisition mode, the triple coincidence be- 
tween two tag scintillators and one Csl detector was required, accompanied 
by at least one good pair of x and y chamber hits. The rare accidental coinci- 
dences involving signals in more than one Csl detector were eliminated from 
the DSF records in the offline analysis. 

The hit anode wire number ua was calculated from 

( Tl -T 2 + N.D) 
ua = 2D ' (1) 



where 7\ and T 2 were the TDC values for two ends of the anode delay line, N 
was the number of wires in a chamber (between 71 and 73), and D was the 
~2.05 ns time delay between the adjacent wires. Truncated wire position xt 
was defined as the nearest integer multiple ua of the 0.8128 cm wire separation 

9- 

x T = Uim(n A -g). (2) 



10 



The anode delay time difference (7\ — T 2 ) depended nonlinearly on the anode 
position. Linearized anode wire positions for each chamber were approximated 
by a polynomial expansion 

m 

x A = Y,c t -(T 1 -T 2 )\ (3) 

i=0 

where the order m of the fit with x 2 P er degree of freedom x 2 /{ m ~ 1)~ 1 
depended on the particular delay line. For our chambers the orders of the fits 
were 3, 4, or 5. The coefficients q were optimized with the program MINUIT [43] 
that minimized the sum of squared deviations (xt — x A ) 2 for each chamber in 
every run. 

The drift distance d x was determined from the drift timing T x using the cali- 
bration lookup table f(T x ) 

d x = f(T x )-T x + AT x , (4) 

where AT X was the drift time offset. The lookup table was calculated assuming 
that incoming cosmic muons are distributed uniformly over the equidistantly 
spaced wires. 

The final hit position x was given by expression: 

x = x T + (-1)™ • si#n(ADC -E + ADC ) • d x + x ff, (5) 

where xoff was the chamber x coordinate offset, and ADCo-e and ADCo 
were the (O— E) ADC signal and associated offset, respectively. 

The individual chamber phases n in the Eq. (5) and absolute horizontal co- 
ordinate offsets xoff and |/off between six chambers were adjusted with 
help of two-dimensional histograms (xi — x 3 ) vs (x 2 — x 3 ) and (yi — y 3 ) vs 
(3/2 - 2/ 3 ) [44,45]. 



6 Raw Data Reduction 

Data summary files were analyzed offline by compressing the DSF event data 
into PAW Ntuples [46]. Drift time-to-distance lookup tables and assorted cham- 
ber calibration constants, as well as Csl detector ADC pedestals and ADC 
temperature corrections, were determined for each run separately. The high 
voltages of the Csl detector phototubes were selected to give ~8 ADC chan- 
nels/MeV scale. The FWHMs of the pedestal peaks were typically 2 ADC 



11 



channels. Therefore, the pedestal widths corresponded to an energy deposi- 
tion of 0.25 MeV, 24 /im pathlength or 0.5% of the most probable energy 
deposition by a minimum ionizing charged particle. Drifts in the pedestal po- 
sition over a one month period amounted to less than 2 ADC channels. 

The precise horizontal coordinate offsets of six Csl crystals inside the dark box 
were determined in the next stage of analysis. Using the preliminary offsets 
read from the plastic template on which the crystals were laying, between 5 
and 10 percent of the reconstructed cosmic muon trajectories were found not 
to intersect any of the predefined detector volumes. After adjusting the x and 
y translation software offset parameters of the Csl modules by maximizing the 
number of tracks intersecting individual crystal volumes with nonzero ADCs, 
the real number of events undergoing scattering in the apparatus or having the 
improperly reconstructed tracks was shown to be below 1 %. That percentage 
was in agreement with the Monte Carlo simulation of cosmic muons interacting 
with the experimental apparatus, showing the most probable muon scattering 
angle of 0.25° and the 2.5° mean scattering angle. The Monte Carlo root- 
mean-square of the scattering angle for the accepted events was 0.66°, Fig. 7. 
That value translates into an average pathlength uncertainty of 1 mm. The 
estimated error in finding the correct crystal position inside the dark box was 
smaller than ~ 0.5 mm, Fig. 8. 

The energy deposited in Csl crystal by cosmic muons along the fixed path- 
length is a broadened distribution due to the statistical nature of energy trans- 
fers. Maximal measured pathlengths were up to 12 cm, while the average path- 
length was ~ 6 cm. A Gaussian model for the energy loss distribution of the 
minimum ionizing particles in Csl thicknesses of less than 12 cm is not a good 
approximation. The distribution of energy losses is asymmetric, with a long 
high energy tail; the most probable energy loss is smaller than mean energy 
loss [47]. 

Light yield temperature coefficients for the individual Csl detectors were deter- 
mined by the "robust" estimation [48] of the ADC values per unit pathlength 
as a linear function of temperature recorded inside the dark box. The least- 
squares condition assumes normally distributed measurements with constant 
standard deviations and is therefore not appropriate in this application, as 
pointed out above. Typically, more than 10 5 cosmic muon events, collected 
over at least one week and spanning the temperature range of ~2°C, were 
used as an input experimental data set. The FORTRAN subroutine MEDFIT, doc- 
umented in the Ref. [49], was adopted as the fitting procedure by imposing 
the requirement of minimum absolute deviation between the measured and 
calculated ADC values per unit path that were dependent linearly on the 
temperature variable. The average light output temperature coefficient for 74 
different Csl detectors extracted by that method was — 1.4± 1.4 %/°C, both for 
fast and total scintillation light components. These numbers are in good agree- 
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merit with the previously reported value of —1.5 %/°C [11], but the spread of 
extracted coefficients for different Csl detectors was large: ±4%/°C (Fig. 9). 
We point out that the method did not allow us to separate the temperature 
dependence of the crystal light outputs from the temperature instabilities of 
the phototubes and high voltage dividers as well as the temperature instabili- 
ties of the ADC modules themselves. The LeCroy catalogue specification [50] 
lists the typical temperature coefficient of a LCR 2249A ADC unit as zero, 
and the maximum coefficient up to ± 3 % for an ADC gate longer than 100 
ns and an average ADC reading of 75 pC (about a quarter of a full 256 pC 
scale) . 

The raw ADC values were corrected for the temporal temperature variations 
and written into the revised DSF files used in subsequent analysis. 



7 Detector Photoelectron Statistics 

The contribution of the photoelectron statistics to the energy resolution of 
Csl crystals was determined by a photodiode-based system. Six individual 
detectors (Csl crystals with photomultiplier tubes and high voltage dividers) 
having identical light emitting diodes (LEDs) coupled to their back sides were 
placed inside the cosmic muon tomography apparatus. The LEDs were pulsed 
at a 10 Hz rate using a multichannel diode driver with continuously adjustable 
output voltage. One split output signal from each channel of the driver gen- 
erated a 100 ns wide ADC gate. The LED light in a Csl detector produced 
fast (~ 20 ns FWHM) PMT anode pulses whose integrated values depended 
on the driving voltage and were equivalent to the fixed cosmic muon energy 
depositions between 10 MeV and 100 MeV. 

A total of five different LED amplitudes were used in measurements of each 
Csl detector. Examples of the LED spectra are shown in Fig. 10. The ADC 
pedestals were recorded simultaneously during the run. The intensity of LED 
light was cross-calibrated against the cosmic muon spectra peaks in Csl crys- 
tals. These muons were tracked in tomography drift chambers and both their 
pathlengths and energy depositions in the crystals could be easily calculated. 
That calculation enabled the establishment of the absolute energy scale in 
MeV. 

The variance a E of a photodiode peak depended upon the mean number of 
photoelectrons N pe on the photocathode created per unit energy deposition in 
the crystal: 

4 = E^ + We, (6) 
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where E was the LED spectrum peak position and of 's were assorted vari- 
ances, such as instabilities of the LED driving voltage and temporal pedestal 
variations. Five measured points (a%,E) were fitted with a linear function (6) 
and the mean number of photoelectrons per MeV N pe for each Csl detector 
was determined. The statistical error of the least-squares fits was typically 
less than one photoelecton/MeV. Using the previously extracted light output 
temperature coefficients all N pc values were scaled to the 18°C point. That 
value was the designed operating temperature of the PIBETA calorimeter. 

The number of extracted photoelectrons per MeV for the fast scintillation 
component fell into the 20-130 range as shown in Fig. 11. The hexagonal 
and pentagonal crystals equipped with three inch phototubes averaged 73 
photoelectrons/MeV, while the half- hexagonal and veto detectors with two 
inch phototubes had a mean of 33 photoelectrons/MeV. The 73/33 ratio is 
very well explained by the (3/2) 2 ratio of photosensitive areas for two different 
photocathode sizes. The measured photoelectron statistics were in agreement 
with the 100 ns ADC gate values of N pe = (20-260) for large (~10cm) pure 
Csl crystals equipped with two or three inch PMTs reported in the past in 
Refs. [15-17,24]. 

The quantum efficiency of our bialkali photocathodes for pure Csl scintillation 
light was 23 % [32] , the average light collection probability for our detector 
shapes 23 % (Sec. 8.2) and the fraction of deposited energy converted into the 
scintillation light was about 12 %. Therefore, starting from the mean number 
of 73 photoelectrons/MeV, we calculated that 1 MeV energy deposited in the 
Csl crystal produced on average about 10 4 scintillation photons. 



8 Monte Carlo Description 

8.1 Tomography System 

Geometrical layout of the experimental apparatus was defined using the GEANT 
detector description and simulation tools [51]. The active elements of the simu- 
lated apparatus were an aluminum box housing six Csl modules, six multiwire 
drift chambers, two scintillator planes and a layer of lead brick shielding, all 
shown in Fig. 4. Generated events were muons with the energy, angular and 
charge distribution of a hard cosmic ray component at sea level. The zenithal 
angle 9 Z distribution of muons at the ground was assumed to be proportional to 
cos 2 8 Z ; the momentum spectrum between 0.1 and 10 3 GeV/c and the energy- 
dependent ratio of number of positive to number of negative cosmic muons 
was parameterized from the data given in Refs. [29,52]. 
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The simulation trigger was defined by requiring a minimum energy deposition 
of more than 0.2 MeV in each scintillator plane. This threshold was about 
one-tenth of the minimum-ionizing peak in the triggered scintillator and cor- 
responded to the discriminator level used in the data acquisition electronics. 

Penetrating cosmic muons and generated secondary particles were tracked 
through the apparatus, and energy depositions, pathlengths in Csl crystals and 
hits in the drift chambers were digitized. All physical processes were turned on 
in the GEANT programs with the cutoff energies of 0.2 MeV for charged particles 
and photons. Inspection of simulated energy deposition spectra showed devia- 
tions from theoretical Vavilov distributions [47] expected in the planar detec- 
tor geometries. The differences were caused by multiple scattering in the ir- 
regularly shaped crystals, and were particularly prominent for the events with 
shorter pathlengths close to the crystal edges. The simulated ADC histograms 
revealed that triggering cosmic muons deposit on average 10.33 MeV/cm in a 
Csl detector, with the most probable energy loss 5.92 MeV/cm. 

Cosmic muon spectra in Csl detectors were described in a satisfactory way 
(x 2 / ( m — 1)~ 1-2) by the combination of a Gaussian distribution and a falling 
exponential tail: 

A(e) = 9(e - e)e-^~^ 2 + 6{e - e )e a - f3e , (7) 



where e is energy (in MeV) deposited in one full Csl module. Parameters e, a e , 
eo, a, and (5 are extracted from the GEANT spectra by imposing the least squares 
constraint to the fit and leaving the pathlength d (cm) as free parameter: 



e{d) = 5.079 + 0.1876d - 9.9390 • lO^d 2 , (8) 

a e (d) = 0.3545 + 7.2329 • lO^d - 3.4148 • 10~ 4 d 2 , (9) 

e (d) = 5.1800 + 0.1989d - 1.0297 • I0~ 2 d 2 , (10) 

a(d) = 4.0581 + 0.4465d - 1.5760 • 10~V, (11) 

(3(d) = 0.7952 + 5.0313 • 10 _1 d - 1.5240 • 10" 3 rf 2 , (12) 



C ^ { 1, otherwise. 

These Monte Carlo spectra, broadened with photoelectron statistics, were used 
to describe our cosmic muon lineshapes produced by the optically uniform Csl 
detectors. 
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8.2 Light Collection Simulation: TkOptics Code 

Propagation of scintillating photons through a uniform detector with ideal re- 
flecting dielectric surfaces and different wrapping materials was studied using 
the TkOptics simulation program [53,54]. The code is a library of FORTRAN 
subroutines with an X- Windows-based user interface written with Tck/Tk 
toolkit [55]. The program can simulate the light output response of an arbitrar- 
ily shaped scintillation detector with given bulk and surface optical properties. 

The polygonal detector shape is defined by the coordinates of its vertices. The 
detector attributes are reflector types of lateral and front crystal surfaces and 
wrapper material, crystal surface-wrapper gap distances and characteristics of 
the photomultiplier tube and a phototube-crystal joint coupling. The program 
handles normal dielectric, specular, diffuse and partially absorbent reflector 
surfaces with arbitrary diffuse fraction, roughness, and specular r s and dif- 
fuse reflectivity r^. Predefined bulk properties of a detector are the medium 
scattering and attenuation length as well as refractive index. A photomulti- 
plier tube is specified by the diameter, position and quantum efficiency of the 
photocathode. 

Different choices of initial scintillating photon distribution are possible, the de- 
fault being a uniform starting distribution throughout the scintillating volume. 
A working volume is a box divided into elementary cubic cells of fixed size. 
Output menu options include the initial and endpoint photon coordinates, and 
direction vectors and timing distributions of the detected scintillation photons 
organized in a PAW Ntuple. Every elementary cell is flagged as a bulk or edge 
crystal cell. Center coordinates of the cells and a fraction of photons gener- 
ated in every cell, and subsequently detected on the PMT sensitive surface, 
are always recorded. Results of high statistics runs with 10 7 photons generated 
uniformly through the detector volume were plotted to show the number of 
photoelectrons as a function of scintillation source position inside the crystal. 
The relative statistical errors of calculated light collection probabilities for the 
bulk crystal cells were < 2 %. 

The bulk attenuation length and the light scattering length of the near- 
ultraviolet emission component inside the Csl crystals were set to 150 cm and 
200 cm, respectively [56] . The index of refraction for the Csl medium increases 
from 1.82 for the blue- green light to about 2.08 for ultraviolet light [56]. 

The simulated light collection probabilities P(x, z) that depend on axial z and 
transverse x positions are shown integrated in the vertical coordinate (y) in 
Figs. 12 and 13. in the form of "lego" plots for one wrapping configuration 
and two different crystal shapes. The coordinate system is defined with the 
z=0 cm origin at the front face of the crystal and a photocathode window at 
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the z=22 cm plane. 

We find that simulated light response of the ideal hexagonal or pentagonal 
PIBETA Csl detector with specular lateral surfaces and a diffuse wrapping 
material can be described with the following parameters: 

(1) The average photon collection probability P for a three inch photocath- 
ode is about 23% for ideal r s =1.0 crystal surfaces and a r d =0.9 diffuse 
wrapper, decreasing by half, to 12 % for r s = 0.8. 

(2) The axial light collection probability variation through the first 10 cen- 
timeters is in the ± 5 % range, with a positive slope dPj dz > for higher 
specular and diffuse reflectivities, namely for r S; d>0.9. 

(3) The axial detected light variation in the back crystal half (z>10cm) is 
up to -30 %. 

(4) The transverse light response referenced to the light yield at the crystal 
axis typically increases towards the crystal surfaces by up to + 5 % for 
z<10 cm, but is generally declining away from the central axis by — 30 % 
at the 2=18 cm plane. 

(5) The root-mean-square of a three-dimensional light nonuniformity is be- 
tween 2.5% and 3.5% for ^<10cm, compared to a ~20% root-mean- 
square for a z> 10 cm volume, where a large spread is caused by the 
inefficient light collection from crystal back corners. 

For lateral surfaces painted with a r^O.9 diffuse substance without an air 
gap we find that the average photon collection probability is lower, ~17%, 
and the axial light collection nonuniformity as a function of axial position is 
always positive and is usually larger than + 10 % in a 22 cm long detector, 
making the root-mean-square of the 3D nonuniformity function > 20 %. 

The Monte Carlo results for the optically uniform half-hexagonal and trapezial 
Csl crystals with the same range of optical properties predict smaller light 
yields and somewhat higher light collection nonuniformities: 

(1) The average photon collection probability with a two inch photocathode 
is about 12 %. 

(2) The typical axial nonuniformity is positive in front, +5%/10 cm, and 
negative in the back crystal half, with a — 30%/10cm variation, 

(3) The simulated transverse light output from the front crystal half is very 
similar to one for the full crystal shapes, increasing up to + 3 % away 
from the detector central axis, but decreasing by almost a third in the 
back corners of the crystal. 

(4) The root-mean-square of detected light output varying between 3.0% 
and 4.5%, in the front (z<10cm) and back (z>10cm) crystal halves, 
respectively. 

These light collection probability distributions calculated for optically homo- 
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geneous crystals explain the major features of the measured optical nonuni- 
formities presented below. 



9 Measured Nonuniformities of Csl Detector Light Responses 

The three-dimensional (3D) spatial distribution of the light output of a scin- 
tillation detector can be specified by giving the number of photoelectrons 
N pe (x, y, z) produced by 1 MeV energy deposition at point (x, y, z) ("3D light 
nonuniformity function"). In the following discussion we limit ourselves to the 
linear one-dimensional variations of the detected light output separable in the 
axial and transverse directions: 



where a z (a t ) is the linear optical nonuniformity coefficient in the z (x) coor- 
dinate, and the coordinate system origin is at the center of the detector front 
face, as explained in Sec. 8.2. 

A simple and straightforward parameterization of the detector light response 
nonuniformity can be made on the basis of scatter plots showing the light 
output per unit pathlength as a function of longitudinal or transverse cos- 
mic muon-CsI crystal intersection coordinates. In this analysis we take only 
cosmic muon events with almost perpendicular trajectories (# 2 >85°), so that 
the measured two-dimensional light outputs are averaged over the pathlengths 
and over ~ 1 cm 2 vertical cross sections. The cosmic muons deposit the en- 
ergy along the well-defined ionization track of the length d in the crystal. This 
description will somewhat underestimate the real light collection probabilities 
(see Sec. 10) due to the integration of a three-dimensional light nonuniformity 
function along the charged particle track. Our work on the fully 3D recon- 
struction of the scintillator light response will be presented in a forthcoming 
publication [57]. 

Fig. 14 shows the axial variation of normalized ADC values as a function of 
distance of energy deposition from the front crystal face for the six represen- 
tative crystals. The axial positions were calculated by averaging two z values 
of the cosmic muon track intersection with the crystal surfaces. 

Panels on Figs. 15 and 16 show the transverse dependence of measured light 
output per unit pathlength, where the independent variable is the distance 
from the crystal axis measured in the horizontal plane. The selected scatter 
plots show the data points for the axial slices at zq=6±1 cm and zo=18±l cm. 




N\ + a z \ ■ z + a t (z) ■ x, z < 10 cm, and x = ±15 cm 
N 2 + a Z 2 ■ z + a t (z) ■ x, z > 10 cm, and x = ±15 cm 



(14) 
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We have chosen to parameterize the axial light output nonuniformity with two 
piecewise linear functions: the light output nonuniformity coefficient a z i(z<10 cm) 
for the front half of the crystal, and the light nonuniformity coefficient a z2 (z>10 cm) 
for the back half of the crystal, both values expressed in %/cm. The transverse 
light output variation at a fixed axial distance zo, a t (z ), is described by the 
average value of a change in luminosity left and right of the crystal axis. 

Several general features are readily noticeable: the spread of measured points 
due to energy deposition straggling, the gradual variation in the collected light 
along and perpendicular to the detector axis, and the decrease in detected light 
when the photon generation occurs close to lateral detector surfaces and, in 
particular, near the back corners of the detector volumes. 

The calculated ADC/pathlength data points shown on the panels of Figs. 14- 
16 have been fitted with linear functions imposing the "robust" condition of 
the minimum absolute deviation between the measured and calculated values, 
Eqs. 14. The normalized light output scatter plots in the axial and trans- 
verse coordinates, as well as two-dimensional ADC/d distributions in the x-z 
bins and associated linear light output nonuniformity coefficients have been 
documented for all studied Csl crystals and are available for inspection at 
the PIBETA WWW site [58]. The average scintillation properties of all the 
measured Csl crystals are listed in Tables 1 and 2. 

In summary, for a set of fifty-nine full hexagonal and pentagonal crystals we 
find: 

(1) The average axial light nonuniformities are a zl {z < 10 cm) =— 0.1 %/cm 
and a z2 (z > 10 cm) = — 1.3 %/cm, respectively (see Fig. 17). 

(2) Distribution of nonuniformity coefficients could be described by a Gaus- 
sian with <r azl ~ 1.3 %/cm. 

(3) The average transverse light output is flat for z<10cm, and the typical 
light variation is —15% at z=18 cm. 

(4) The Kharkov-grown crystals on average have twice the axial optical nonuni- 
formity of the Bicron-grown crystals. 

For fifteen half-hexagonal and trapezial crystal shapes we find: 

(1) The average axial nonuniformities are a z i(z<10 cm) = — 0.3 %/cm and 
a z2 (z>10 cm) = -1.0 %/cm. 

(2) The transverse light output variation is limited in the front crystal half, 
but increases up to — 30 % in the back of the crystal. 

With the insights gained in the simulation of the light collection probabili- 
ties (Sec. 8.2) we conclude that the ±2 %/cm spread of the nonuniformity 
coefficients corresponds to an equivalent range of 0.8-1.0 in crystal surface 
reflectivities, and/or 100-250 cm range in Csl attenuation lengths. 
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10 Light output uniformity scans of Csl crystals with 137 Cs gamma 
source 

A light-tight plywood box was made to house a single Csl detector and as- 
sociated measurement apparatus described below. A 0.662 MeV 137 Cs gamma 
source was placed on a moveable support next to the detector and collimated 
by a 20x10x5 cm 3 lead brick with a 6 mm collimator hole. The photomul- 
tiplier analog signal from the detector was processed with an ORTEC 454 
timing amplifier with a gain of 30 and a 50 ns integration time constant. The 
amplifier output was digitized with a peak sensing ORTEC AD811 ADC unit. 
The same signal was discriminated and produced the trigger rate of about 
5 kHz with a 137 Cs source present. The background rate, without the source 
present, was 50-100 Hz. 

Pedestal runs were taken with a clock trigger and used to properly correct 
offsets of the ADC spectra. Temperature variation during one run was less 
than 0.3°C, typically causing the overall light output variation of less than 
0.5 % in a ADC spectrum that was gated with a 100 ns window. Several 
runs were taken with the source removed to find the shape of the background 
spectrum. It was confirmed that the background spectrum does not depend 
on the position of the lead collimator, so the same background lineshape was 
used in the analysis for all source positions. 

The collimated 137 Cs source was placed by remote control in turn at five points 
along the axis of each crystal at 2, 6, 10, 15 and 20 cm from the front face of 
the crystal. 

Following background subtraction, recorded spectra were fitted with a sum 
of a Gaussian and exponential function, Fig. 18. The peak position and the 
FWHM for each spectrum were extracted with the statistical uncertainty lower 
than 0.3 %. 

The dependence of the peak position on the placement of the source is illus- 
trated in Fig. 19 and follows the trend of the tomography results. The GEANT 
simulations revealed that the 662 keV gamma rays could probe all the volume 
of Csl crystal, but the shower energy deposited at the central axis is only about 
one-sixth of the energy deposition near the crystal surface. The cosmic muons 
transfer the energy uniformly along their tracks in Csl material: the scintil- 
lation volume over which the measured light output is integrated is therefore 
larger, and inherent averaging leads to smaller extracted optical nonuniformity 
coefficients. This feature is borne out in the panels of Fig. 18 where the to- 
mography data for three different crystals yield consistently smaller axial light 
nonuniformities when compared with the radioactive source measurements. 

The simple 137 Cs scans with the described apparatus were used to evaluate 



20 



the light collection nonuniformities for up to six crystals per day. 

Radioactive source scan measurements similar to ours are described in Refs. [22,31]. 



11 A GEANT simulation of the PIBETA calorimeter response to 10 
120 MeV positrons and photons 

We have studied the simulated PIBETA calorimeter response to the 69.8 MeV 
-K + ^e + u positrons and to the (69.8+2m e +) MeV photons, where m e + is the 
positron rest mass. The goal was to find the intrinsic difference between the 
responses to monoenergetic positrons and photons with identical total incident 
energies, emphasizing the correct modeling of the low-energy tail below the 
edge of the Michel [/,— >e + vv spectrum at 52.8 MeV. 

The 0.46 % accuracy of the new TRIUMF measurement [4] of the ir + — >e + v 
branching ratio was limited by event statistics and systematic uncertain- 
ties in evaluating the low-energy tail of the positron peak. The 7r + —>e + v'-f 
positrons and the forward-peaked bremsstrahlung gammas were detected in a 
46-cm-diameterx51-cm-long Nal(Tl) crystal "TINA". The tail correction of 
1.44± 0.24 % for the energy region from MeV to 52.3 MeV was determined by 
subtracting the measured Michel positron spectrum from the detector positron 
response functions measured with the 20-85 MeV e + monoenergetic beams. 
An additional tail component of ~ 0.4 % due to radiative processes was es- 
timated by Monte Carlo method. The TRIUMF group apparently made no 
attempt to account for the potential difference in the scintillator light response 
of positrons and photons and neglected potential nonlinearities of the energy 
scale. Both effects could arise from the light collection nonuniformities of the 
Nal(Tl) detector. 

In the recent PSI 7r + ^e + z/ experiment [3] the quoted 0.29 % systematic uncer- 
tainty of the extracted branching ratio was also dominated by the 1.64± 0.09 % 
electromagnetic loss fraction as well as by the 0.95± 0.19 % contribution from 
the photonuclear reactions. The ~4-7r calorimeter used in that measurement 
was made of 132 hexagonally shaped BGO crystals, with the light yield claimed 
to be homogeneous to within 1.5 % over the whole 20 cm length of every crys- 
tal. 

The lineshapes of positrons and photons in our PIBETA calorimeter and their 
dependence on the optical properties of 240 constituent Csl crystals were pre- 
dicted using the GEANT detector description and simulation tools. The code 
defined the geometry and tracking media for the complete PIBETA detec- 
tor [59] . The inner detector region was occupied by two cylindrical wire cham- 
bers and a segmented plastic veto detector. The incident particles, monoener- 
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getic positrons and photons, were generated in the center of the crystal sphere. 

Photonuclear reactions in the active detectors as well as in the surrounding 
passive material were modeled in a user subroutine [60] that was added to the 
default GEANT version 3.21. The probabilities of photonuclear reactions were 
calculated using the published cross sections from the reaction thresholds 
up to the energy of 120 MeV [61-66]. 

The 240 individual Csl module shapes were specified, taking into account 
the mechanical tolerances of the physical Csl crystals, with Teflon and alu- 
minized Mylar wrappings filling the 200 /xm intermodule gaps. The irregular 
Csl modules had to be constructed from up to six GEANT generalized trape- 
zoidal wedges. Every Csl crystal volume was considered a sensitive detector 
with the associated luminosity N pe expressed as the number of photoelectrons 
per MeV, and two axial and one transverse light collection nonuniformity coef- 
ficients, a z i, a Z 2, and a t , respectively. A set of seventy-four detector luminosity 
and nonuniformity values, extracted in the tomography analysis, were initial- 
ized in the GEANT-accessible database. The optical parameters of the remaining 
163 modules were drawn randomly from the distributions in Figs. 11 and 17. 
The assumption was that the crystals produced in the future will be of the 
same optical quality as the ones that are already delivered. 

We have always simulated at least 10 5 events for every fixed set of the calorim- 
eter parameters (N pe ,a z i,a Z 2,a t ). The low-threshold trigger was defined by re- 
quiring a sum of the ADC readings exceed 5 MeV of the light-equivalent energy 
for one calorimeter supercluster containing 32 crystals. The integrated detec- 
tor acceptances with a low-level trigger were 85.8% for ~ 70 MeV positrons 
and 85.5% for ~ 70 MeV photons, so all extracted quantities had the rela- 
tive statistical uncertainties of ~0.2%. The response of the calorimeter was 
parameterized by the FWHM of the simulated ADC spectrum and the tail 
contribution being between the adjustable low level (default LT being 5 MeV) 
and high level thresholds (two default values HTi and HT 2 being 54 MeV and 
55 MeV, respectively). 

For the purposes of comparison and gain normalization we first studied the 
response of an ideal, homogeneous calorimeter, with realistic Csl crystal light 
outputs. The "software" gain factor of every individual Csl detector was de- 
termined from the fitted peak positions of the simulated ADC spectrum sum 
over the crystal with maximum energy deposition and its nearest neighbors. 
The peak value positions were found after smoothing the histograms with a 
multiquadric function, to improve on the limited simulation statistics. The 
ratios of peak positions in the ideal, homogeneous calorimeter and the corre- 
sponding values for the physical, nonuniform calorimeter, were, by definition, 
the individual detector gain corrections. The values of these ratios were refined 
in three steps of the iterative procedure. The extracted energy resolution of 
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the nonuniform detector in principle depends on the convergence of the gain- 
matching process and the resulting uncertainties of the individual detector 
gains. That gain-matching process was equivalent to the PMT high voltage 
adjustments in the real experiment that are effected to obtain the best en- 
ergy resolution with the modular detector. We estimate that our simulation 
procedure fixes the software detector gains with the accuracy of ~ 1 %. 

The set of gain constants depended on the optical properties of our Csl crystals 
as defined in our GEANT database, but also upon the incident particle chosen 
for the Monte Carlo calibration runs and its energy because of differences in 
shower developments of photons and positrons. The gain normalizations of 
the calorimeter detectors were calculated by aligning the simulated -K + —>e + v 
positron peak positions. The same procedure for matching the Csl detector 
gains was used in the real data-taking runs. Due to the gain renormaliza- 
tion, the influence of transverse light nonuniformities smaller than 10 % on 
the ADC spectrum peak positions and the ADC lineshape could be neglected 
in comparison with the axial light nonuniformity effects. The simulated ADC 
values were summed over all detectors with the energy deposition above the 
1 MeV TDC threshold. We also examined the calculated ADC sums for the 
clusters which contained the crystal with the maximum energy deposition and 
its nearest neighbors. The full-widths at half maximum for these spectra are 
labeled in the following figures and tables as FWHM( 2 2o) an d FWHM(nn) , re- 
spectively. No event-to-event uncertainties of the ADC pedestal values were 
assumed in the simulation. Stability of the PIBETA electronics tested under 
real experimental conditions in 1996 calibration runs and the quality of the al- 
gorithms used for the first and second pedestal correction reduces the pedestal 
peak root-mean-square to ~8 channels, equivalent to ~ 0.3 MeV. 

The GEANT-calculated resolution of the PIBETA calorimeter consisting of 220 
Csl detector crystals and 20 veto crystals was parameterized by the fractional 
full width at half maximum: 



FWHM (220 ) /rwN 2.36 • <7 e 



(%) = -^r^ = (15) 



E e + E, 



(54.69 + 1.197£ e+ + 0.4656 • 10~ 2 E 2 e+ ) / E°J mo ^ 



\ (37.19 + 0.060£ e+ + 0.2734 • 10~ 2 E 2 e+ ) / E°+ 5256 
for the positrons with the most probable peak energy E e + MeV and 



FWHM (220) _ f (37.07 + 0.5507£ 7 + 0.3233 • 10~ 2 E 2 ) / £°- 6549 
E y ^ ' ~ \ (40.52 + 0.9437 • £ 7 + 0.4931 • 10" 2 £ 2 )/£ 7 



(17) 



for the monoenergetic photons with the detected peak energy E 1 MeV (Fig. 20). 
The top lines refer to the cases of an optically uniform detector (a z =0) while 
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the bottom equations describe the predicted response of the nonuniform cal- 
orimeter (a z =T0M0GRAPHY) . 

The average fractional energy resolution FWHM( 2 2o) of 5.3 % (6.0 %) was 
achieved for the ~ 70 MeV kinetic energy positrons in the optically uniform 
(nonuniform) calorimeter, as compared to 5.7% (6.8%) resolution for the 
equivalent energy photons. That was the resolution found with no cuts ap- 
plied on the light-equivalent energy generated in the calorimeter vetoes. Re- 
quiring less than 5 MeV detected in the veto shield decreased the statistics 
by about 10 % and suppressed the low-energy tail, but did not improve the 
energy resolution at the peak position. 

Limiting the ADC sums to the nearest-neighbor crystals (six or seven crystal 
clusters) the fractional FWHM( NN ) energy resolution could be parameterized 
by: 



FWHM (NNW a _ / (69.19 + 3M7E e+ + 0.4198 • 10~ 2 E 2 e+ ) / E°J 560 
E e+ [ } ~ \ (68.13 + 2.479£ e+ + 0.3368 ■ 10~ 2 E 2 + ) / E°J 018 

for the positrons with the detected peak energy E e + MeV and 



FWHM(nn) _ / (43.16 + 2.318£ 7 + 0.4148 • 10~ 2 E 2 ) / E°- 7723 
" E 7 1 j ~ \ (39.57 + 2.302£ 7 + 0.1329 ■ 10~ 2 E 2 ) / E°- n " 

for the monoenergetic photons with the detected peak energy E 1 MeV (Fig. 21). 

The average FWHM( NN ) for the 70 MeV positrons and gammas in nearest- 
neighbor nonuniform crystal clusters was 8.7 % (6.1 MeV) and 9.9 % (6.8 MeV), 
respectively. Imposing the 5 MeV cut on the detector veto signals improves 
these resolutions only marginally. These numbers should be compared with the 
response width of the ideal uniform nearest-neighbor clusters: 8.2 % (5.7 MeV) 
and 8.7% (6.2 MeV). 

The percentage of the events in the tail between the preset low and high energy 
threshold was tracked in the same GEANT simulation. 

The listing of the low-energy tail contributions for 69.8 MeV e + and 70.8 MeV 7 
in Table 3 shows that in the optically uniform calorimeter they differ by 1.8 %. 
The optical nonuniformity does not change the positron tail contribution, but 
increases the photon tail by ~ 0.2 %. The imposition of the 5 MeV veto-shield 
cut decreases the electromagnetic leakage to 2.4% for positrons and 4.7% for 
photons, a e + -7 tail difference of +2.3%. 

Table 4 shows the low-energy tail components for the simulated nearest- 
neighbor cluster ADC sums: with no applied cuts the contributions were 4.7% 
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and 6.7 % for e + 's and 7's, respectively. The veto shielding cuts decreased both 
fractions by about 0.3 %. 

If the positron originates from a radiative ii + ^e + h"y decay, its average Monte 
Carlo peak position is essentially unchanged. Its tail contribution is unaffected 
if simulated ADC sums extend over full calorimeter but are increased by ~ 
1.4% if only the nearest-neighbor clusters are summed. The radiative decay 
matrix elements used in the calculation were taken from Ref. [67]. All quoted 
numbers have an approximate systematic uncertainty of ~0.2%. 

The detector response to positrons and photons in the energy range between 
lOMeV and 120 MeV is also nonlinear because the volume distribution of 
the shower energy deposition depends on the incident particle type and the 
energy. Our GEANT studies showed that the magnitude and energy dependence 
of these nonlinearities does not change significantly because of Csl crystal 
light collection nonuniformity. Results of the calculations, both for the case 
of a homogeneous and optically nonhomogeneous calorimeter, are displayed 
in Figs. 22 and 23. Nonlinearities of the detector response in the covered 
energy range for positrons and photons are very close in magnitude and shape 
of energy variation and amounted to ~1.5% for the ADC sums over 220 
crystals and up to ~2.8% if the simulated ADC sums were restricted to the 
over-the-threshold ADCs of nearest neighbors. 



12 Conclusion 

We have measured the optical properties of seventy-four pure cesium iodide 
crystals that were polished and wrapped in the diffuse Teflon reflector. The 
results are summarized separately for the full-sized and half-sized crystals in 
Tables 1 and 2. 

The deduced light yields parameterized by two axial and one transverse light 
collection nonuniformity coefficients constitute a minimum set of parameters 
necessary for a realistic Monte Carlo simulation of the modular Csl calorime- 
ter. 

The predicted energy resolutions FWHM( 2 2o) f° r ~ 70 MeV positrons and pho- 
tons in the full PIBETA calorimeter with the ideal, optically uniform Csl mod- 
ules with specified luminosities were shown to be close, 3.7 MeV and 4.0 MeV, 
respectively. The upper limit of the low energy tail contributions in the re- 
gion between 5 MeV and 55 MeV were calculated to be 6.9% and 8.7% for 
the positrons and gammas in an optically homogeneous detector, respectively. 
After applying the 5 MeV calorimeter veto cut, these tail corrections decrease 
to 2.2 % and 4.5 %. 
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The average deduced axial light nonuniformities of real Csl crystals wrapped 
in a Teflon sheet had negative slopes, — 0.18%/cm and — 1.6%/cm, for the 
front and back half of the crystal volume, respectively. The corresponding 
GEANT simulation of the nonuniform PIBETA apparatus for 68.9 MeV e + and 
70.8 MeV 7's shows that both energy responses will be broadened to the av- 
erage FWHM(22o) of 5.9 % and 6.9 % while the associated tail contributions 
will change only for photon spectra, increasing the low-energy tail by ~0.3%. 
The simulated calorimeter ADC spectra are shown in Figs. 24 and 25. 

The nonlinearity of the measured energy scale caused by the optical nonuni- 
formity is <2.8% throughout the relevant e + /7 energy range. This spread is 
consistent with the precision of energy calibration required to extract the tail 
corrections with the systematic uncertainty of ~0.2%. 

The predicted ADC spectra of the monoenergetic positrons, electrons and 
tagged photons in the energy range 10-70 MeV will be compared with the 
measured responses of the partial Csl calorimeter arrays in a forthcoming 
publication [68]. 
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Fig. 1. The PIBETA Csl shower calorimeter consists of 240 tapered hexagonal, 
pentagonal and trapezial pyramids. The crystal ball has an opening for beam entry 
and a symmetric one for light readout from inner detectors (not shown). 

Fig. 2. The dimensioned drawing of the HEX-A (top) and HEX-D1 (bottom) Csl 
crystal shapes. 

Fig. 3. The fast-to-total light component ratio for the all studied Csl crystals. The 
light outputs from optically polished and unwrapped crystals were measured imme- 
diately after the crystals were received from a manufacturer with a digital oscillo- 
scope at the average ambient temperature of 22° C. 

Fig. 4. The Monte Carlo simulation of penetrating cosmic muons interacting with 
the experimental apparatus. Six Csl crystals were defined inside a dark box and 
deposited energies and track lengths of cosmic muons and generated secondary 
particles in detectors were digitized. 

Fig. 5. The schematic diagram of the electronics logic for the tomography apparatus. 
The calibration event was triggered by two-fold scintillator coincidence SlxS2. The 
trigger used during cosmic muon data acquisition required three-fold coincidence 
between two scintillators and at least one Csl detector, (SlxS2)xJ] Cslj. 

Fig. 6. One-dimensional residuals defined as differences between the best straight 
line fits from the reconstructed coordinates in the hit chambers 1 and 2 and the 
measured positions in the top chamber 3. The FWHM of the peak is 0.9 mm for 
x coordinate and 1.5 mm for y coordinate, respectively. The 1.2 cm wide plateaus 
contains ~ 5 % events with incorrect left-right identification and /or the events un- 
dergoing the large angle scattering in the apparatus. 

Fig. 7. The Monte Carlo distribution of scattering angles of the cosmic muons 
intersecting a HEX-A Csl crystal shape. The root-mean-square of the scattering 
angle for the events with the reconstructed tracks is 0.66°. That value corresponds 
to the Csl pathlength uncertainty of ~ 1 mm. 

Fig. 8. The average measured cosmic muon pathlength in Csl crystals was ~6cm. 
An accurate determination of the track length requires the precise knowledge of 
the absolute crystal position with respect to the drift chambers. The number of 
muon tracks intersecting a Csl crystal volume is shown as a function of horizon- 
tal crystal offsets. Contours of the 25-level plot reveal a crystal position with an 
uncertainty of ±0.05 cm (left panel, pentagonal Csl crystal: £ ff = 12.32± 0.05 cm, 
z oS = 30.50± 0.05 cm). The same type of plot is shown for one hexagonal Csl crystal 
on the right panel. 
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Fig. 9. Fast light output temperature coefficients for all analyzed Csl calorimeter 
detectors. The mean value is — 1.4%/°C, with the temperature coefficients ranging 
from — 4%/°C to +4%/°C. Temperature changes inside the light-tight box (filled 
markers) and average ADC values (open circles) for one detector with the negative 
temperature coefficient of — 1.9/°C are shown in the inset. Data points shown in 
the inset span one week of data acquisition. 

Fig. 10. LED spectra in one Csl detector recorded for five different values of the LED 
driving voltage. The horizontal scale is a raw ADC count of the fast Csl scintillation 
component (a 100 ns ADC gate) corrected for the pedestal offset. The equivalent 
energy range covered is 10-100 MeV. The photoelectron statistics can be determined 
by measuring the width of a photodiode peak a\ as a function of a spectrum mean 
value E. The linear dependence described by Eq. (6) is demonstrated in the inset. 
The inverse slope of the straight line equals the number of photoelectrons created 
per unit energy deposition. 

Fig. 11. The number of photoelectrons per MeV of deposited energy for all studied 
Csl calorimeter modules. The left panel shows the distribution for the full-sized 
hexagonal and pentagonal Csl detectors viewed with three inch phototubes. The 
right panel is the equivalent histogram for the smaller half-hexagonal and trapezial 
detectors equipped with two inch phototubes. 

Fig. 12. The simulated light collection probability as a function of axial (z) and 
transverse (x) coordinates for an optically ideal hexagonal detector (HEX-A) with 
a two-layer Teflon wrapping. The size of the histogrammed two-dimensional bin is 
lxl cm 2 . The average light collection probability with a three inch phototube is 
23%. 

Fig. 13. The Monte Carlo light collection probability as a two dimensional function 
of axial and transverse coordinates for an ideal half- hexagonal detector (HEX-D1) 
with a Teflon wrapping. The average light collection probability with a two inch 
photocathode is 12 %. 

Fig. 14. Relative light per unit pathlength ADC/ti as a function of axial position 
of energy deposition is shown for six Csl crystals. The phototubes were always 
positioned on the surface with the larger values of the z coordinate. Plotted values 
were corrected for temperature variation effects and the PMT gain drifts. 

Fig. 15. Detected light per unit pathlength as a function of transverse coordinates 
of energy deposition shown at the fixed distance along the crystal axis (zo=6± 1 cm) 
for six Csl detectors. 

Fig. 16. Detected light per unit pathlength as a function of transverse position 
of energy deposition shown at the fixed axial distance (zo=18±lcm) for six Csl 
crystals. 
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Fig. 17. The light nonuniformity coefficients a z \ and a z2 (%/cm), determined sep- 
arately for front (z <10cm) and back (z >10cm) crystal sections, are shown here 
in a scatter-plot for 74 Csl detectors. The full-sized hexagonal and pentagonal de- 
tectors are represented with the full marker points, while the open circles indicate 
the half-hexagonal and trapezial detector shapes. The dotted line is the result of a 
"robust" straight-line fit: a z i(%/cm, z< 10 cm)=-l. 44+1. 19-a z2 (%/cm, z>W cm). 



Fig. 18. The 662 keV 137 Cs gamma source and the background spectrum measured 
at the source axial position z=10cm from the front face of the pentagonal S003 
crystal. 



Fig. 19. The comparison of normalized light output (ADC reading/pathlength) 
along the axis of three different Csl crystals deduced from the cosmic muon to- 
mography data (full line), and as measured with the 137 Cs gamma source scan 
(dotted line). 



Fig. 20. The predicted fractional energy resolution FWHM( 220 ) /E(220) °f the PI- 
BETA calorimeter for positrons and photons in the energy range 10-120 MeV. The 
GEANT-deduced relative FWHM (220 ) for the sum over 220 Csl ADCs is shown both 
for an ideal homogeneous detector (full lines) as well as for a nonuniform calorimeter 
(dotted lines). Notice that the zero value on the vertical scale is suppressed. The 
curves represent the parameterizations quoted in the text. 



Fig. 21. The GEANT-deduced FWHMjnnj/^nn) f° r the sum over nearest-neighbor 
Csl ADCs is shown both for an ideal homogeneous detector as well as for a nonuni- 
form calorimeter. 



Fig. 22. Nonlinearity of the PIBETA calorimeter energy response for positrons and 
photons in the energy range 10-120 MeV. The ratio of the detected ADC peak posi- 
tion to the incident particle energy is calculated in a GEANT simulation. Results for 
an ideal homogeneous detector and a nonuniform calorimeter are shown separately. 
The normalization is given by ADC/£'(70MeV e + ,j)=l. 



Fig. 23. Nonlinearity of the PIBETA calorimeter energy response for positrons and 
photons in the energy range 10-120 MeV. Ratio of the detected peak position of 
summed nearest-neighbor ADC spectra to the incident particle energy is calculated 
in a GEANT simulation. 
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Fig. 24. The predicted PIBETA calorimeter spectra of the monoenergetic 69.8 MeV 
positrons and 70.8 MeV photons. The 220 ADC values above the 1.0 MeV TDC 
threshold were summed. A GEANT simulation assumed different numbers of photo- 
electrons/MeV and linear axial light collection nonuniformities as extracted in the 
tomography analysis. The results for the "realistic" detector depend on the treat- 
ment of the Csl crystal surfaces and the details of the gain matching algorithm. 

Fig. 25. The predicted PIBETA calorimeter spectra of the monoenergetic 69.8 MeV 
positrons and 70.8 MeV photons. The ADC value for the Csl crystal with the max- 
imum ADC value and the ADC readings of its nearest neighbors were summed. A 
GEANT simulation assumed different numbers of photoelectrons/MeV and linear ax- 
ial light collection nonuniformities for each detector, as measured in the tomography 
apparatus. 
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Table 1 

Average scintillation properties of hexagonal and pentagonal PIBETA Csl calorim- 
eter shapes PENTAs, HEX-As, HEX-Bs, HEX-Cs and HEX-Ds (59 crystals). All 
crystals were polished and wrapped in two layers of a Teflon membrane and one 
layer of aluminized Mylar film. The light yields are normalized to the tempera- 
ture of 18°C. All other parameters were measured at the average laboratory room 
temperature of 22° C. 

Bicron Csl Kharkov Csl 
Crystals (22) Crystals (37) 



Fast-to-Total Ratio (100 ns/1 fjs gate) 0.835 

# Photoelectrons/MeV (100 ns ADC gate) 83.4 

# Photoelectrons/MeV (1 ^s ADC gate) 99.9 
Fast Light Temp. Coefficient (%/°C) -1.20 
Total Light Temp. Coefficient (%/°C) -1.25 
Axial Nonuniformity Coefficient (%/cm), 

z<10cm, 100 ns ADC gate -0.200 
Axial Nonuniformity Coefficient (%/cm), 

z>10cm, 100 ns ADC gate -1.62 
Transverse Nonuniformity Coefficient(%/cm), 

z =6cm, 100 ns ADC gate -0.400 
Transverse Nonuniformity Coefficient(%/cm), 

z =18cm, 100 ns ADC gate -1.10 



0.739 
65.1 
88.0 
-1.61 
-1.40 

-0.124 

-1.61 

-0.500 

-1.20 
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Table 2 

Average scintillation properties of half-hexagonal and trapezial PIBETA Csl calo- 
rimeter shapes HEX-Dl/2s and VETO-l/2s (15 crystals). All crystals were polished 
and wrapped in two layers of Teflon foil and one layer of aluminized Mylar. 

Bicron Csl Kharkov Csl 
Crystals (3) Crystals (12) 



Fast-to- Total Ratio (100 ns/1 fis gate) 

# Photoelectrons/MeV (100 ns ADC gate) 

# Photoelectrons/MeV (1 ^s ADC gate) 
Fast Light Temp. Coefficient (%/°C) 
Total Light Temp. Coefficient (%/°C) 
Axial Nonuniformity Coefficient (%/cm), 

z<10cm, 100 ns ADC gate 
Axial Nonuniformity Coefficient (%/cm), 

z>10cm, 100 ns ADC gate 
Transverse Nonuniformity Coefficient(%/cm), 

zo=6 cm, 100 ns ADC gate 
Transverse Nonuniformity Coefficient(%/cm), 

z =18cm, 100 ns ADC gate 



0.806 
34.5 
42.8 
-2.08 
-1.78 

-0.326 

-1.00 

-0.755 

-4.10 



0.728 
31.5 
43.3 
-0.63 
-0.69 

-0.500 

-1.83 

-0.827 

-5.14 
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Table 3 

The predicted energy resolutions and tail contributions for 69.8 MeV e + and 
70.8 MeV 7's events in the full PIBETA calorimeter. The light output in photo- 
electrons/MeV and the linear axial light collection nonuniformities measured for 
individual Csl crystals (a,zi );Z 2=T0M0GRAPHY) were used in a GEANT simulation. The 
values for the perfect optically homogeneous crystals (a z i, 2 2=0) were shown for 
comparison. 



Parameter 


69.8 MeV e+ 




70.8 MeV 7 




No Ey cut 


dzl,z2=0 


TOMOGRAPHY a zl)Z2 


a z i, Z 2=0 


TOMOGRAPHY a z i >z2 


Peak Position (MeV) 


68.87±0.03 


68.90±0.03 


70.10±0.03 


69.97±0.03 


FWHM (220 ) (MeV) 


3.66±0.03 


4.10±0.03 


3.98±0.03 


4.76±0.03 


5<Events<54MeV (%) 


6.44±0.09 


6.46±0.09 


8.20±0.10 


8.47±0.10 


5<Events<55MeV (%) 


6.89±0.09 


6.98±0.09 


8.86±0.10 


9.11±0.11 


E v <5MeV 


a z i,z2=0 


TOMOGRAPHY a zl)Z2 


a z i,z2=0 


TOMOGRAPHY a zl>z2 


Peak Position (MeV) 


68.88±0.03 


68.92±0.03 


70.07±0.03 


69.99±0.03 


FWHM (2 2o) (MeV) 


3.67±0.03 


4.09±0.03 


3.99±0.03 


4.76±0.03 


5<Events<54MeV (%) 


1.91±0.05 


2.02±0.05 


3.92±0.07 


4.15±0.07 


5<Events<55MeV (%) 


2.24±0.05 


2.36±0.05 


4.46±0.07 


4.68±0.07 
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Table 4 

The predicted energy resolutions and tail contributions for 69.8 MeV e + and 
70.8 MeV 7's events in the PIBETA clusters containing the crystal with maximum 
energy deposition and its nearest neighbors. 



Parameter 


69.8 MeV e+ 




70.8 MeV 7 




No E v cut 


a z i,z2=0 


TOMOGRAPHY a zl , z2 


«zl,z2=0 


TOMOGRAPHY a zl , z2 


Peak Position (MeV) 


67.74±0.03 


67.85±0.03 


68.80±0.03 


68.65±0.03 


FWHM (NN) (MeV) 


5.22±0.03 


5.90±0.03 


6.01±0.03 


6.80±0.03 


5<Events<54MeV (%) 


3.83±0.07 


4.01±0.07 


5.47±0.08 


5.78±0.08 


5<Events<55MeV (%) 


4.54±0.07 


4.73±0.08 


6.33±0.09 


6.67±0.09 


E v <hMe\ 


a z i,z2=0 


TOMOGRAPHY a zl , z2 


a>zi,z2=0 


TOMOGRAPHY a zl , z2 


Peak Position (MeV) 


67.75±0.03 


67.83±0.03 


68.80±0.04 


68.64±0.03 


FWHM (NN) (MeV) 


5.52±0.03 


5.83±0.03 


5.92±0.03 


6.76±0.03 


5<Events<54MeV (%) 


3.53±0.07 


3.73±0.07 


5.28±0.08 


5.59±0.08 


5<Events<55MeV (%) 


4.19±0.07 


4.39±0.07 


6.10±0.09 


6.42±0.09 
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